Soliton structures in a molecular chain model with saturation 
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00 \ In the present work, we study, by means of a one-dimensional lattice model, the collective excita- 

, tions corresponding to intra molecular ones of a chain like proteins. It is shown that such excitations 

' are described by the Nonlinear Schrodinger equation with saturation. The solutions obtained here 
are the bell solitons, bubbles, kinks and crowdons. Since they belong to different sectors on the 

j^jrj, parametric space, the bubble condensation could give place to some important changes of face in 

r- { ' this kind of nonlinear system. Additionally, it is shown that the limiting velocity of the solitons is 

' the velocity of sound waves corresponding to longitudinal vibrations of molecules. 

\Q . PACS numbers: 05.45-a,87.10.+e 

(N 

Ph ' I. INTRODUCTION 

d 

In the present days the research to the dynamic properties of one dimensional molecular chains have increased. 
The structure of a great amount of macromolecules is represented by subunits with mutually weak, slightly flexible 
bounds connecting them with each other. The important biological structures of this type are RNA, proteins and 
: DNA polymer chains. The peculiarity of bio-polymers is that, they are heterogeneous, and their elementary subunits 
^ ' have complex structures and carry long-lived nonlinear excitations. As it is well known the propagation of energy 
and electrons in protein molecules is the crucial factor for maintaining the life of biological systems. So, the problem 
of storage and transportation of energy through protein chains arises. The energy used in biological cell comes from 
\ the energy of liberation during the process of hydrolysis of adenosinetriphosphatc (ATP) molecular structures. The 
. energy of this process is of approximately 0,42 eV (or 2250 cm~^). Proteins consist of chains of hydrogen-bonded 
\l ' peptide groups, three of these chains in a helical arrangement define the a— helix structure In his seminal work, 
Davydov proposed an explanation of the fundamental transportation problem of energy released by hydrolysis of 
' adenosinetriphosphatc and transferred to proteins in biological systems ^] . This energy remains localized and moves 
. . , along the protein chains at a reasonable rate to perform useful biological functions. It could be trapped and transported 
^ ' in proteins as quanta of the intra molecular C=0 stretching mode, the so called amide-I vibration, with excitation 
energy around 1650cm^^. The localized spatial region where the energy is trapped can propagate along the protein 
chain, in such a way that a soliton-like mechanism for energy transport is possible. This problem of transporting 
d I energy from one point to another inside the cell is a long-standing problem that remains of great interest. 

Besides, from experimental point of view, we can also discover a lot of contributions that are directly related to a 
similar phenomenon in the DNA. For instance, in the experimental work on short DNA rings e.g., Refs. 0, II the 
kink tendency of DNA sequences were studied. In the last years the great amount of works devoted to the nonlinear 
dynamics of DNA shows that this area is an active field, for example theoretical proposals of wrapping DNA around 
the nucleosome, where kinks play a great rule, were proposed in various papers, see the Refs. [^,[^,13]. Concerning the 
a— helical protein dynamics, some works have been dedicated to study this system including high order excitations 
and different molecular interactions, in the discrete and continuum level, [1, 

In this contribution we investigate soliton-like structures within the framework of a certain generalization of the 
Davydov's model, considering the case of neighboring interactions as the same class before and after a peptide group. 
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In the next section we briefly expose the modified Davydov's model leading to the Hamiltonian that will be used in the 
later. Section III is devoted to derive the nonlinear cubic - quintic Schrodinger equation by a suitable transformation 
from the original nonlinear Schrodinger equation with saturation. The soliton structures of this equation with some 
specific characteristics are presented in section IV. Finally, in the last section we deliver some comments. 



II. DAVYDOV'S MODEL. 



Due to its transparency and seminal properties, Davydov's model continues to encourage intense work regarding 
the research of the nonlinear treatment of molecular systems 2] . In his pioneer works he and co-workers demonstrated 
that the corresponding nonlinear equations for the molecular excitation in the quantum treatment admit solitonic 
structures. Indeed they assumed that the energy transportation in proteins is carried out by means of transportation 
of amide - 1 vibration. In what follows we consider an infinite chain of weakly bound molecules (or groups) with a mass 
M and a distance R from each other. Internal excitations of molecules (electronic or vibrational) are characterized by 
an energy and an electric dipole moment d directed along the chain. The internal excitation of molecules and their 
motion around equilibrium positions are inseparably linked. 

In the case of a one-dimensional chain, only interactions between neighbor molecules are taken into consideration. 
Below, we follow the modeling done by Davydov and his co-workers. If the intramolecular excitation has the energy 
e, the collective excitations of this model can be described by the Davydov's Hamiltonian: 



n 

+T + U, (1) 

where index n labels the molecule that occupies position r„ in the chain, while y _B„ are creation and annihilation 
of intramolecular excitation boson operators. The quantity J — 2cPR^^ characterizes the transition of intramolecular 
excitation due to the resonant interactions while d is the electric dipolar moment. The last two terms in ([!]), as usual, 
correspond to the kinetic and potential energies of the longitudinal displacements. 

When the i-th molecule is excited, the static interaction with neighboring molecules of this molecule changes. This 
is reflected by the introduction of the function _D„. The displacement p„ from the equilibrium distance R in the state 
I ^) is deflned by the expression 



p„ i?- (r„ - r„_i). (2) 

In the state | 0) without intramolecular excitation, the chain has periodicity and the intermolecular distances are R. 
Let us now consider the function D„, which in the nearest neighbors interaction limit has the following form 



Dn = Sri (k«-l - fn|) + (kn - T-K+ll) 

where D = 2J)„ (i?) , /?, 7 are parameters of the theory. 

The potential energy of the molecules in the non excited chain is chosen in a harmonic approximation under the 
assumption that the constant spring uj is the same for all of them. In this case we can express the potential energy as 



n 

and the kinetic energy can be written as 



n n \ — oo>/<n / 
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where dot over the letter represents temporal derivative, p = According with the quantum mechanics treatment 
the collective interactions of interest can be described by the wave function 



n 

where coefBcients ipn {t) are normalized as ^ \ tpn (t)\^ = 1- These coefficients characterize the distribution of exci- 

n 

tations along the molecular chain, tpn being the probability of finding the quantum system in the site n or \tpn 
being the density of probability of finding the excitation. The equation for determining these wave functions can be 
obtained from the Schrodinger equation 



zft^|*)=H|*), 

that can be reduced, using the exphcit form of the operator H, Eq. ([T]), the fact that the functions |0) correspond 
to different values of n and are orthogonal each other, to obtain the following system of equations 



ih 



d_tpn 

dt 



e + T + U 



n+l J ■ 



Pi 



D 



(4) 



The functional i.e. the hamiltonian that can be associated with this equation of motion can be written as F 



e+T+U- 1 



R 2i?2 



D 



Following Toda [10|, it is convenient to associate the displacements pn with their canonically conjugate variables 

l>n 



Sn — — —M ri. Then the kinetic energy can be expressed in terms of these new variables as 



T 



2M 



In the next section we will see the derivation for the nonlinear Schrodinger equation (NSE) with a saturation term. 



III. NSE EQUATION WITH SATURATION 



Now, we can derive the equation of motion for the displacements and for their canonical conjugate variables s„, for 
this we consider that 



9F 1 

asn M 

opn n 

After eliminating the variables s„ from the preceding system of equations, we find the equation for the displacement 



{2pn - Pn+l - Pn-l) 



M 



1^: 



n+l| 



\i'n-l 



(^^" '^"'^ ~ IV'n+ll^ - Pn-l I'f/'n-ll^) 



(6) 
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The system of equations ( |4]) and ([6]) defines the collective excitations and deformation of a chain. 

For an analytical treatment if we are interested on the distribution of excitations along the chain, we turn to the 
analysis in the continuum limit. For doing this, let us introduce the dimensionless variable £, = and the continuous 
fimctions as usual p{^,t) and {^,t) such that 

p (n, t) = pn (t) , ip {n, t) = ipn (t) ■ 

Expanding p ± 1) and ± I) in series in the standard manner 



Pi^±i,t) « pio± 



2 de ' 



d 



1 



\rPi^±i,tr « i^(^,t)r±^i^(^,<)r + -^iv4e,i)r, 

and keeping only terms up to second order of magnitude, we transform the equations ([4]) and ([6]) to the following 
system of two equations: 



and 



ih 



dt 



d^p{^,t) ^d^p I3D 



+ 



p^D ( a2|^|2 dpd\i^\^ , ,,2' 



M 



+ =0, 



(7) 



with 



A 

T + U 



e + T + U -D~2J 



I 00 / 00 



(8) 



and Va = (w/M)2 , where VaR = Va is the acoustic longitudinal velocity in the chain. 

We will look for traveling solutions moving along the chain with some velocity V — vR. In this case, we have the 
following transformation 



P (e, t)^p{^^vt), ^ (e, i) = - vt) exp {id (e, t)} . 
Replacing the equation ^ in equation ([7]) and after integrating, we obtain 



(9) 



P - 



(3D 



I3D I ; |2 I ' 



(10) 



with G= {vl- v"^) . 

Finally, we substitute the equation (|10p into the equation (jIII|) and obtain the nonlinear equation 



ih 



dip 



dt 



J 



ki ivi 



k2 iv^r 



V' = o, 
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with the values fci = ^j^tjj, ^2 = 1 2R*m'^ ' ^3 = mr'^ • Rewriting these parameters in terms of the exciton-phonon 
coupling constant x — we have ki = ^ , k2 = 2RM^ • Equation (IIII|) is the well known NSE with saturable 
nonlinearity. This equation arose earlier in various branches of physics, particularly in nonlinear optics and simulated 
saturation (decrease) effects of the nonlinear response of a medium in large electromagnetic fields [l3l | 

Let us further simplify equation (jIII[) . If we take into account in this equation, the nonlinearity not higher than 
0(6|?A|^), and fc2 >> ^1^3, we obtain for the distribution of excitations the Cubic - Quintic Nonlinear Schrodinger 
Equation(CQNSE) 

As known, nonlinear equations similar to equation (fTTjl possess interesting structures when the attractive and 
repulsive terms could compensate each other. So, in the next section we report some solutions that appear as natural 
excitations along the molecular protein chain. 



IV. SOLITON STRUCTURES 



For solving the equation of motion presented in the previous section we have to consider physically boundary 
conditions. Since we are interested on the fact that the displacements of the perturbed units could only take the local 
character, it is proposed that at long distances from the occurring perturbations, displacements are very weak and 
practically the distribution of excitations at long distance vanishes i.e at " infinity" it is zero. The second boundary 
condition is considered when the displacements take constant values at infinity. These restrictions of our chain at 
"infinities" could be fixed for the time evolution of the perturbations along the chain. 



A. Trivial boundary condition. 

For simplicity, let us analyze the case when fc2 < and ki > 0. If this is done, the last term in (|lip represents the 
repulsive part of the nonlinearity and the 4th term the attractive one. Further, if we make the variable transformation 
= i ' ^ = yf^' = = ~X2MRG^ ■^^^ finally obtain: 

^^ + ^-,^+{m'--m')i^ = 0. (12) 

The CQNSE ([T21) was studied from various points of view, here we follow the results and conclusion obtained in 
the works [lli,[3and 

The corresponding solution of the CQNSE with trivial boundary condition 

ip for X —>■ ±oo, 

is the so called drop-type soliton that is not a topological soliton because the vacuum also has the same asymptotic 
value. This implies that for the equation we have the static non-topological soliton [ll[. The non topological 
soliton solutions are those of which boundary conditions at infinite are the same that vacuum state. However, 
topological solitons have boundary conditions different from the vacuum. This means in particular, that states of 
degenerated vacua might exist. It is to say that soliton "will be moored" by its boundary conditions. An example of 
topological soliton solution is a step or kink. For our case we have the drop soliton 



16a 



ij} = e*'''V-4Q: ( 1 + + — cosh {^/~^{x - xq)) 



(13) 



with 



a — —^v 



A7 

2(3D' 



(14) 



The traveling soliton should be obtained using the Galileo transformation in the same form what we define a traveling 
wave in mechanics 



FIG. 1: Displacement of the molecules represented by the equation ( llSp around the equilibrium position for A — 2.5, a = 
0.1,7 = 2.3, 1.2 



V v'^ 

cosh [\/—a {x — a;o)) cosh [\/—a {x — Vt — xo)) . 

The soliton solution (jl3p has the normalized motion integral named the "number of particles" JdxjV'l' = 1, 
calculating this integral we find 

16a _ 1 ^ 
~ ^ cosh[l/V3] 

The approximate value for the parameter a ~ —0.051. Replacing this value in the relation p4p we obtain the 
restrictions of the main parameters ^ = 0.102. The quintic part of the nonlinear equation produces the effect of 
counterbalance the attractive forces between "two particles " in the mechanical analogy method represented by the 
cubic nonlinearity. 

If we replace (fT3| in (fTO]) we have the distribution of changes in the relative distance between molecules. 



p{^, t) = (^-^] 7- . r^^,, (15) 

V <J-fR J l + v cosh {£. - ^o)\ 



with rr - 1 + 4«B7P _ 1 2A7= _ 

Wltll cr - 1 + jy/G - 1 - -MG^ ^ - a ' 

Some numerical representations of this solutions with different values of the main parameters are represented in 
the figure [2 

In the particular case when the parameters satisfy the relation 

1 At' 



3PD MG V MG 
we have the solution 



'^^''^"i?fl'^2A£)cosh^[iV^(e-Co)]' ^''^ 



which is the new distribution of changes in the relative distance between particles. 
From (jl5ll6p the maximum deviations di and d2 are 
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From these equation we can see that the displacement due to the appearance of sohtons corresponding to Eq. ([15]) is 
greater compared with the similar equation PB]) . This means that the strong " damage" will be caused by the soliton 
represented by the equation (fT6|l . We can see, that the presence of solitons leads to the pronounced deviation of the 
peptide groups from their equilibrium position and the second one should produce a breakage of the chain system 

when 2^ = 1. 

Using the expression ([5]), we can obtain the values of the total energy of the peptide group displacement as follows: 




B. Condensate Boundary conditions 



We can suggest also, that there are very specific restrictions that could cause soliton excitations to appear along 
the protein chain. Besides the natural or well known bell soliton excitation it is also very possible the appearance of 
other types of solutions. For example, there could surge topological or non topological solitons because the CQNSE 
supports them. It is well known that the equation (|lip supports kinks and bubble type of solitons. For this case, we 
will use the well reported results obtained by many authors specifically we mention [11]. 

Let us see the case of regular solutions of the CQNSE pi|l in 1+1 dimensions of the space-time. We rewrite this 
equation in a slightly different form by using the ground states and putting them in the equations. In order to visualize 
the ground states, we will use the following form of the CQNSE: 



+ 'f,, - (3 - {2A + ao)) {\^\^ - do) = 0. (17) 

This version permits us to find the soliton solutions in explicit form. It can be demonstrated that the eq. (|17p 
could be generated from the relation with the help of the following scale transformations 



r 



Making some algebra, parameters A and cto are related each other by the relation 

:+^-—(i + ^/T^^). (19) 

4 /it/ v / 



A „ 3 1 

— = -2-- 

(To 4 /it 



Without loss of generality it is possible to fix the value of Uo = 1, because properties of the solutions depend easily 
on the parametric relation {A/uo). Here the parameter A can be both positive or negative based on the physical 
requirement we could impose on the system. Further, the variables C and r will be treated as if x and t were the 
customary variables. 

1. Bubble solitons. 

For obtaining gray or bubble solitons when the degenerated vacuum is not absolute, it is imposed the boundary 
condition 

X ±00 and if K3. 



with K3 — \/ . For this case, the solution of equation (fTTll takes the form 
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FIG. 2: This displacement is similar to the well known "classical bubble" solution in field theory, but for our case it represents 
the agglomeration of units in the molecular chain and for this reason it is named "crowdon" ( I2ip 



eWcosh[2V^(^-yo)] 



with y = C, — VT and 



This bubble is a nontopological solution, and their topological charge is equal to zero, Qh = 0. Here the parameter A 
satisfies 1 < A < A. then for the displacement p according to the equations (fTU|) and (|20p we obtain the expression 



27 



cosh^ [2Va(y - yo] 



ai + a2 sinh [2^{y - yo] 



(21) 



being 

{A -I) 



and 



while 



"i="^ + 9^(2^+l)(4-A)' 



3G 3G 
4- A ' 2A+1 



.= I^,„..^^^(^. (22) 

In the case when the ai.2 are both positive, it is easy to see that ai > cti since for this case 4 > ^ > 1. The 
displacement corresponding to the bubble like soliton excitation Lp\, is a typical gray soliton and can be depicted in 
Fig [21 Apparently this type of solution is similar to others obtained for nonlinear classical models. But in contrast to 
the well known feature, in our case we have not a bubble displacement, instead we have an agglomeration of molecules 
that travels along the chain i.e. we have here typical crowdon solution forming the agglomeration of molecules. This 
agglomeration is traveling along the chain like an accumulation of molecules conserving velocity and profile. 

This result is linked to solutions that are moving slowly with less velocity in comparison with the sound velocity, 
i.e. when v^ — v^ > 0. But when the opposite occurs, i.e. when the soliton velocity is greater than the sound velocity 
we have soliton solution on the step or on the background. This type of solution is presented in the Fig. [31 



FIG. 3: Typical soliton like perturbation on the background that represents the displacements. The interpretation that could 
be given to this solution is that in a real manner represents a bubble or the region inside of which the displacements are 
increasing and this effect is contrary to the agglomeration. 



So, the sign of G determines the type of soHton solution for the displacement that corresponds to the distribution 
of excitations. When 1 < yl < 4, bubble solitons for the excitations appear corresponding to displacements of both 
types: popular gray solitons and normal "bell" soliton on the background. In the first case for displacements we have 
obtained the crowdon solution. In the last case, this displacement shows a huge separation between neighborhoods 
in the molecular chain. This solution travels along the protein with velocity v greater than the acoustic velocity. 
Of course, we have here a little difference of popular bubbles of other nonlinear equations. As it is well known, the 
classical bubble is characterized by the rarefaction inside crucial region, but in this case, the feature of our solution 
is contrary to this best known property of classical bubbles. Indeed, the displacements outside the perturbed region 
have constant values, while inside of this region the separation between molecules increase as shown in Fig. [3l When 
this displacement is strong enough, one could observe the breaking of the protein chain and the solution could be 
transformed to a peak or cusp like non-classical solution. 



2. Crowdon solutions 



The "space" of our model is a line with two points as boundaries at the infinite. The existence of topological soliton 
named kink, is due to the properties of the mapping between the degenerated minimums of the potential, with the 
space, saying better, with its boundaries, that in this case it is a discrete set. 

Kink solitons appear when the potential supports degenerated absolute vacua in = K3 = ^^^^ , when A > 4. So, 
the boundary conditions for this case, are 

(C = -00) = -K3, ip{C = oo) = K3. 
The kink solution has the following form 

fk = J= , (,2dj 



V3 Jl + ^cosh2[2V^(y-j;o)]' 



For antinkink, as usual, we can take the inverse signs, so when x —t +00, the field will approach the value — K3,. 
Calculating the topological charge of this solution we obtain — 1- For the antikink, we have Qant = — 1- In this 
case, the regular solutions with finite energy are divided in 4 topological sectors. The sectors with finite energy can be 
characterized by means of the following pairs of indices (— K3, K3), (ks, — K3), (— K3, — K3) and (^3, K3) that correspond 
to the values of the field at infinities i.e. for <p (a; = —00) , [x = -l-oo). The kink, antikink and trivial solutions 
Lp (x) = ±K3 are the members of the 4 sectors. 

Using the equations ((23|) and (fTO| For the displacement p{C,, r) we have the following expression 

^ m(3D / sinh^ [^Va{y - yo)] \ ^24) 
' RM U3+«4Cosh2[2V^(2;-yo)]y ' ^ ' 



10 




FIG. 4: The evolution of the displacements shows typical characteristic of "crowd" soliton solution. For a qualitative repre- 
sentation of the solution, we take the values 04 = 1,03 = 3,u = —1.2. Apparently, this picture represents "classical bubble" 
solitons but inside the perturbed region the displacements reduce their values forming the crowdons (|24p . 



with 



"3 



3G 



A 



3G 



2A+1 



and K is defined in equation ((22l) . 

Considering positive values for when the boundaries values of displacements are fixed, there is a possibility of 
emergence of crowdons that will propagate along the chain. The displacements show the typical picture of " classical" 
bubble solitons. This bubble in the molecular chains context is not a dip in the background. On the contrary, this 
solutions represent the wave agglomeration of molecular units traveling along the chain. So, along the molecular chain 
the distribution of excitations evolves like a kink soliton while its corresponding displacement evolves like typical crowd 
like soliton solutions, see Fig. [H 



C. Stability and Self- localizations 

Next, we can use some results concerning the possibility of stabilizations of those types of solitons we have found 
in the section A and B for those distribution of excitations with non- topological charges. We will follow the findings 
of the paper [l^ 

As it is well known, our system is closely related to the phi-six theory, that shows several important behavior in 
nuclear hydrodynamics, ferromagnetism, phase transitions and other branches of physics and natural sciences [llj. 
The equation of motion that we will analyze in the context of self localization is the equation \12\ obtained above. 

The main ideas of the approach developed in 11 1 and (isj could be summarized as follows: This method in some 
sense is similar to Lyapunov functional approach for analyzing stability problems. For the general multidimensional 
case of spherical or cylindrical symmetry the dynamical system associated with the solitons NSE is: 

r d<P 

where $ = iIj*'4' , the dimension of the space is V and with suitable change of parameters we use the potential part 
of the energy as [ll|, [lB| 



U ($) ^ [/($) = $ + 2A$ - (2 + .4) $2 + $3 



with 



A ^ -2 - —{1 + VT+4a), 
4a 
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and a is defined by the equation p4|) . For our case, a molecular system, the equation ((25|) is the same equation (fT2|) 
with V = \ [l2| • The constant solutions of the equation Ijlip at A = undergoes a supercritical bifurcation. 

For analyzing the self-localizations of soliton structures we will evaluate the behavior of the following functional 

B{t) 



B{t) = tt2 



■D-l 



(26) 



where N denotes the "number of particles" integral of motions. The magnitude B{t) represents in some sense the 
width of the soliton solution obtained in this model. As it can be seen from the equation ([26|) the functional B{t) is 
proportional to the expectation value of the square of the position. From here this magnitude is proportional to the 
width of the soliton solution. Let us evaluate the second derivative of this magnitude and we obtain 



d^Bjt) 



Z TT 



8 / \tPr\ dr + AV 



— $ -U]dr 



(27) 



As we are concerned with solitons in one dimensional case for the molecular system, finally making V = 1, and 
evaluating the conditions under self-localization of soliton structures we finally obtain that following conditions 



liV=l, ^<0, for A<2, (28) 

where 2J denotes the "volume" of the soliton structure. Now let us calculate the possible values of the parameters for 
getting self localizations of solitons. We make a = —g and for satisfying the condition (pSj) we obtain the inequality 

i (l + v/r^) < y (29) 

with g = As we can see from the relation for all values of g < 0, the self-localization is possible i.e. when 

the parameters satisfy this relation < 0. In the case of positive region of the axis p, we can conjecture that 
self-localizations would occur when 

That means, we should try to find stable soliton like solutions within this sector that could determine the possible 
values of the relevant parameters in real experiments. 



V. CONCLUSIONS 



We have taken into consideration in some sense an improved Davydov's model, and obtained a nonlinear evolution 
equation for the displacement that in mutual coordination with the excitation generate nontopological and topological 
solitons. 

When the molecular system evolves, it is of course interesting to know which types of structures it could support 
during its evolution when its dynamic is modeled by the cubic quintic nonlinear Schrodinger equation. When the 
boundary condition is of the trivial one, we obtained a typical well known "bell" soliton that represents the separation 
of the molecular units inside the perturbations. Let us make some comments about the physical meaning of the solution 
obtained when we considered the condensate type of boundary conditions. In this specific case, the molecular model 
subjected to both "condensate boundary conditions" , in both directions, the displacements p acquires a constant value. 
This means that due to multiple internal or external factors that affects the molecular system, the displacement far 
from the perturbed region acquires some certain values that in some sense could be considered as established during the 
process of the existence. In the case, when bubble like solutions appear as natural excitations, displacements around 
the central part of the perturbed region could be of two types: the first one corresponds to the crowd solutions, that 
means we have in some sense an agglomeration of units that compound the lattice i.e. a "crowdon" . This agglomeration 
travels along the chain with velocity Vg less than that of the sound. The second solution for displacements is the 
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"soliton on the step" and travels faster than sound. It represents the increasing of distance between elemental units 
inside the perturbed region of the molecular chain. For the case, when the excitations reveal the typical profile of kink 
solutions, as we can see from equations (fM)) . the displacements apparently evolve like "classical bubble" solitons but 
they represent the crowdons again, and represent the agglomeration of units along the molecular chain. According 
to the stability analysis, we also conjecture that it is possible to obtain stable solitons inside the narrow segment in 
the parametric domain determined by the equation (jIV C|) . Finally, it should be interesting to check in more realistic 
models, the behaviors of these structures after the inclusion of , for example, dissipation term to describe the effects 
from water molecules surrounding the molecular system. 
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